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Using extreme-ultraviolet attosecond-pulse-trains, we investigate the photoionization dynamics 
of a Helium atom in the presence of moderately-strong lO^^Wcm"^) femtosecond laser pulses. 
The electronic structure of a laser-dressed atom is traced in real-time through precision measure- 
ments of ion-yields and photo-electron angular distributions. Quantum interferences between photo- 
excitation paths are interpreted using the Floquet formalism. As the laser pulse intensity ramps on 
femtosecond timescales, we observe transitions between ionization channels mediated by different 
atomic resonances. The quantum phase of interfering paths is extracted for each channel and com- 
pared with simulations. Our results elucidate photoionization mechanisms in strong-fields and open 
the doors for photo- absorption/ionization control schemes. 



The recent advances in 'attosecond science' have given 
a new impetus to the study of atomic and molecular phe- 
nomena by providing direct real-time access to electron 
dynamics [1. Experiments in this regime are typically 
conducted using extreme-ultraviolet (XUV) attosecond 
pulses or pulse trains along with precisely synchronized 
strong- field femtosecond near-infrared (IR) laser pulses, 
to obtain new insights into dynamics of electronically ex- 
cited systems[lH3]. As the roots of 'attosecond science' lie 
in the strong-field concepts developed in 1990's [4 , the 
application of new attosecond techniques to refine our 
understanding of atomic/molecular dynamics in strong 
fields is of particular interest [1. 

Here, we report precision real-time measurements of 
the transient non-equilibrium electronic structure of He- 
lium in intense fields. We investigate the quantum in- 
terferences in two-color photo-ionization pathways us- 
ing XUV attosecond pulse trains (APT) and variable 
strength near-infrared (IR) laser fields. As the field in- 
tensity changes on femtosecond timescales, we observe 
switching between ionization channels. We find that 
yield from each resonance-mediated ionization channel 
oscillates with a specific phase. We interpret this quan- 
tum phase using Floquet interaction model. Numerical 
calculations using time-dependent Schrodinger equation 
(TDSE) serve to elucidate the important role of Floquet 
interferences in photo-excitation and ionization. 

We use amplified 65fs, 785nm IR pulses of 1.5mJ en- 
ergy, which are split into two parts. One part is focused 
onto a Xe gas-filled hollow waveguide to generate APT 
of high-harmonics (HH). The APT along with the co- 
propagating driving infrared pulse (IR^) is focused onto 
a He gas jet using a toroidal mirror. The second part, a 
probe pulse (IRp), goes to a delay stage and is focused 
on the He target with a 50cm lens. The schematic of our 
experimental set-up is provided in [5 . Photo-electrons 
are imaged using a velocity-map-imaging (VMI) setup. 
He+ ions are spatially imaged such that there is a one- 
to-one correspondence between the point of origin and 
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FIG. 1. (a) XUV spectrum and the relevant He states (b) 
Calculated XUV photoabsorption cross-section as a function 
of laser intensity, (c) Floquet manifold showing one-photon 
spaced components of IR-dressed 2p state and the two inter- 
fering paths associated with HH13 and HH15 excitation. 



the point they hit the detector, allowing us to eliminate 
Gouy phase averaging and obtain a high-quality signal [5]. 

Fig. [TJa) schematically shows the experimental HH 
spectra relative the unperturbed He resonances. The 
15th harmonic is resonant with the 5p electronic state 
and 13th harmonic is slightly below 2p resonance. Two 
other harmonics that we observe, i.e. 11th and 17th, are 
much weaker (20 times lower) and non-resonant, hence 
they do not play a significant role in this study. In Fig. 
[TJb) we show the photo-absorption cross-section of He as 
a function of photon energy and peak IR intensity calcu- 
lated using the method described in [6 . Clearly, the dis- 
crete atomic resonances evolve into a complicated struc- 
ture even in a moderately intense laser field of the order of 
lO^^Wcm"^. The higher excited-states (3p, 4p, 5p, etc.) 
exhibit positive shifts, which can be approximated by the 
ponderomotive effect. The low- lying 2p state exhibits a 
negative shift and develops multiple branches. At inten- 
sities around 5xlO^^Wcm~^, the excited state structure 
bears very little resemblance to the unperturbed case. 
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In previous experiments with APT and IR, it has been 
observed that the He+ ion yield oscihates as a function 
of time delay with half-IR-cycle periodicity [JHS]. This 
has been interpreted as interference between wave pack- 
ets generated by successive bursts in the APTjT) Hm [TT]. 
Using Floquet interaction picture [12 , it has been shown 
that the interference between different Fourier compo- 
nents of an IR dressed state leads to this oscillatory 
variation of ionization probability at 2uj frequency [9l [13]. 
Importantly, recent theoretical work [13 has raised cru- 
cial questions about the role of different ionization paths 
and the phase of ionization signal that have not been 
addressed by the experimental studies conducted so far. 
Here we measure the phase of the ion-yield oscillations 
and investigate its variation in terms of the evolution 
of laser-dressed atomic structure. We demonstrate that 
for a given IR-dressed atomic resonance (Fig. [ijc)), the 
quantum phase difference between the dipole transitions 
to different Fourier components determines the oscilla- 
tion phase of the ionization signal. 

In our XUV+IR ionization measurements, the IR field 
results from a combination of two IR pulses. The weaker 
pulse iJRd) is phase- locked to the APT and the stronger 
pulse (IRp) is time delayed relative to the APT[5 . For 
two main harmonics (HH13 and IIII15 in our case), the 
probability of ionization by the XUV and combined IR 
field can be obtained for a given Floquet state as [13 
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where Mq{A{t)) and M2(A(r)) are the IR-intensity de- 
pendent matrix elements representing transitions to the 
direct and two-photon-dressed Fourier component. The 
terms /o and /2 are the strengths of two harmonics, uj 
is the central frequency of the IR, r is the time delay 
between XUV and probe IR, is the phase at which 
the attosecond pulse is locked to the driver IR field. The 
quantum phase difference between Mq and M2, which is 
an important quantity in the paper, is represented by 
(j). Figure [TJc) diagramatically shows the Floquet mani- 
fold associated with 2p resonance and the two interfering 
excitation paths discussed above. 

The delay dependence of combined IR ampli- 
tude and phase, i.e. A{t) and ^(r), is given as 
A{t) = {Al^ Al^2ApAdCos{ujT)f^^ and 3{t) = 
sm~^{Adsm{ujr)/A{T)) , where Ap and Ad are the am- 
plitudes of the probe and driver IR fields. For a weak 
driver-IR field {Ad <C Ap), the dominant frequencies in 
P(r) are luj and 2uj. We define the normalized amplitude 
of ion-yield oscillation as Posc{r) = {P{r) - Pavg)/Pavg, 
where Pavg is the one-cycle average. The Pqsc in the weak 
driver case can then be approximated as 

Posc(r) = Ti cos(cjr) + T2 cos(2cjr + 25q + 0) (2) 

where Ti , T2 are the amplitude factors. 
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FIG. 2. (a) Normalized He^ ion-yield oscillations for XUV 
and 2IR pulses (driver and probe) at different probe peak 
intensities (o.c. denotes optical cycles). Asymmetry in the 
double peak oscillation structure is reversed between 1.4 
TWcm~^ and 3.4 TWcm~^. (b) Fourier amplitude of ion- 
yield oscillations shows two prominent frequencies. 



Figure |2|a) shows normalized experimental IIe+ ion- 
yield as a function of the time delay between APT+IR^^ 
and probe IRp at different probe IR intensities. The 
driver IR intensity is less than lO^^Wcm"^. We observe 
a distinct oscillation structure with one-cycle and half- 
cycle components in accordance with equation (2). The 
oscillation at lu frequency, which arises from the inten- 
sity modulation due to the IR-IR interference, acts as a 
reference with respect to which we can robustly measure 
the phase (j) of the Floquet path interferences occurring at 
2uj. The 2IR method thus allows use to compensate for 
the interferometric drifts and other experimental varia- 
tions. Furthermore, this approach also provides a method 
to measure So, which represents the timing of attosecond 
bursts relative to the peak of driving IR field. 

Fig. |2|b) shows the Fourier transform of ion-yield os- 
cillations over the time-delay range that spans more than 
20 optical cycles. At low intensities, the loo component 
dominates, however, as the probe intensity is increased, 
the 2uj component increases. This is expected as the 
IR amplitude modulation decreases and the two-photon 
dressed Floquet contribution (i.e. M2 in equation (1)) 
becomes increasingly important with intensity, leading 
to a stronger interference signal. Importantly, the os- 
cillatory structure in Fig. [2|a) at two probe intensities, 
namely, 1.4 TW cm~^ and 3.4 TW cm~^ is very differ- 
ent. The asymmetric double-peak structure at 1.4 TW 
cm~^ shows the left-peak to be higher, whereas, at 3.4 
TW cm~^ the situation is reversed and the right-peak 
is higher. This difference in oscillation structure is a di- 
rect manifestation of the change in phase relationship be- 
tween the 2uj and luj components. The intensity depen- 
dent change in relative phase originates from the change 
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FIG. 3. (a) Experimental XUV+IR photo-electron spectrum 
of He at different IR probe intensities. As the IR inten- 
sity is increased the 2p + 3cj peak becomes dominant, (b) 
Calculated photo-electron spectrum (c) Experimental VMI of 
photo-electrons at different IR probe intensities. 



in value of (j) in eq.(l), which represents the phase differ- 
ence between the two interfering contributions. Before 
we quantitatively discuss this phase, it is important to 
identify the Floquet states contributing to the ionization. 

To identify the Floquet paths, we utilize photoelectron 
spectroscopy. Fig. |3ja) shows the experimental electron 
spectra at probe intensities used in Fig. [2] The observed 
electron peaks in Fig. [sj^a) are associated with IR ioniza- 
tion of XUV excited 5p, 4p and 2p atomic states. At low 
intensities, the ionization is mediated by 5p resonance, 
and we observe a strong peak corresponding to 5p + luj 
process. This is expected as 15th harmonic is initially 
resonant with 5p state. As the intensity is increased, 
the Ap + luo starts contributing. At higher intensities the 
2p+3co' channel dominates the ionization signal. This ob- 
servation is in accord with Fig. [ijb) as the 2p structure 
Stark shifts downward in energy with increasing inten- 
sity, becoming resonant at higher intensities. 

TDSE calculations for comparable intensity parame- 
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ters also yield similar results (Fig. 13 
resolved photoelectron images in Fig. 3 
as intensity is increased towards 3.4 TWcm~^, side lobes 
corresponding to the 'g-wave' structure appear. This an- 
other indication of 3-photon ionization of XUV excited 
state 2p state. Thus, the results from Fig. |3] confirm 
that the dominant two-color ionization pathway changes 
from the 5p-mediated ionization at low intensities to 2p- 
mediated ionization signal at higher intensities. Next, we 
extract the phases for various resonance mediated ioniza- 
tion channels and establish quantitative relationship be- 
tween the strong-field variation of atomic structure and 
intensity dependent phase change observed in ion-yield 
oscillations of figure |2j 

Figs. [2] and [3] show that even moderately intense laser 
pulses can significantly modify an atom and its ioniza- 
tion dynamics. Now we probe these transient dynamics 



in real time. Fig. [Zj^a) shows the raw He+ ion yield 
as a function of time-delay at 3.4 TW cm~^ probe in- 
tensity. The negative time-delay axis implies that the IR 
probe arrives ahead of the XUV pulse. As the time-delay 
changes from -20 o.c. towards zero, the asymmetric 2uj 
oscillation structure develops and evolves from high-left- 
peak asymmetry to high-right-peak asymmetry. Note 
that this behavior is similar to Fig. [2j except that Fig. 
|2] shows the dependence of ionization signal on peak in- 
tensity, whereas. Fig. [4]^ a) explores the dependence on 
intensity variation within the IR pulse profile. 

Next, we extract the phase of 2uj component of ion- 
yield oscillations relative to the luo component using stan- 
dard Fourier-transform methods. Figure [Ij^b) plots this 
phase as a function of the time-delay (solid line). We 
repeat the phase extraction excercise on the ionization 
signals in Fig. |2|a) and obtain the phase of 2uj compo- 
nent at different peak intensities. Assuming a Gaussian 
profile, we calibrate the time-delay axis such that each 
delay value corresponds to a specific instantaneous inten- 
sity (top- axis of Fig. pFb)). Thus, we can plot the phases 
extracted from Fig. ^2|a) as square dots in Fig. |4]^b) 
and compare them with delay-dependent phase curve. A 
good match between the two independent results demon- 
strates the soundness of our experimental method. 

The insets in figure gb) show the TDSE results for 
XUV+IR ionization oscillation at 3.4 TWcm~^ and 1.5 
TWcm~^ with arrows pointing to comparable experi- 
mental points. At 3.4 TW cm~^, where 2p contribution 
strongly dominates, the phase of 2uj oscillations is zero 
(right-inset), implying ionization yield maximizes when 
APT arrives at the peak of IR- field. We can use this 
fact to extract the 5q in eq. (2) and remove APT tim- 
ing offset from our phase deduction, thereby setting the 
starting point of experimental phase value to zero. With 
this adjustment, the modified form of eq. (1) becomes 
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Starting from zero, as we move to longer time-delays 
(i.e. lower intensities), the experimental phase value in- 
creases. The phase goes through tt and reaches 1.57r at 1.4 
TWcm-2 intensity. The TDSE calculation (left-inset) 
also yields a similar result as the 2uj oscillations in shifted 
in phase by ~ 1.57r at 1.5 TWcm~^. 

The variation of phase figure [4]^b) can now be under- 
stood in terms of the quantum mechanical phases asso- 
ciated with various ionization channels and the relative 
dominance of different channels. The calculated energy- 
resolved photoelectron data in figure [4]^c) elucidates the 
role of different channels. At high intensity (3.4 TW 
cm~^), the 2p resonance mediated ionization channel at 
0.69 eV dominates. The oscillation of this channel has 
zero phase and the ionization yield peaks at zero-delay 
where XUV pulses arrive at the peak of IR field. For this 
to happen the quantum phase (j) between the Mq and M2 
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FIG. 4. He^ ion yield as a function of time-delay at 3.4 
TWcm~^ probe intensity, (b) The phase of the 2u component 
of ion-yield oscillations with time delay, (c) TDSE results 
showing the energy-resolved oscillations in electron yield of 
various resonance mediated ionization channels. 



turn phases associated with XUV+IR ionization chan- 
nels, which are not known a priori. These observations 
are general and should be valid whenever XUV excitation 
occurs in a strong field. As the modification of electronic 
states depends on the instantaneous value of the IR laser 
intensity at the time of XUV excitation, tuning IR in- 
tensity can change contribution of the different matrix 
elements to the interfering terms. By precisely control- 
ling the XUV spectrum, IR intensity and the time-delay 
between the APT and IR, it is possible to control the ion- 
ization dynamics. A recent experiment |14j demonstrat- 
ing XUV transparency validates some of these ideas. 

In conclusion, the two-color photo-ionization of He by 
XUV and IR pulses represents a relatively simple yet rich 
system for exploration of light-matter interaction. We 
show that even a moderately strong field, often used in 
XUV-IR pump-probe studies, can substantially modify 
electronic structure and ionization dynamics on attosec- 
ond timescales. We identify different ionization channels 
and measure the quantum phases associated with these 
Floquet paths, which is the first such measurement to our 
knowledge. Understanding the evolution of atomic struc- 
ture in strong fields and its interaction with attosecond 
XUV pulses can provide us with a knob for fine control 
of photo-dynamics. 

This work was supported by the National Science 
Foundation (NSF) under contract PHY-0955274. 



transition matrix elements in eq. (3) has to be zero. In 
other words, 13HH and 15HH induced transitions to the 
Floquet components associated with 2p excitation are 
'in-phase' (Fig. [ijc)). 

As intensity decreases, the 4p channel at 0.72 eV starts 
contributing substantially. Interestingly, the oscillations 
in this channel are completely out of phase with the 2p 
contribution (Fig. [4|c),1.5 TW cm~^ ). Using eq. (3), 
this implies that = tt for 4p mediated ionization chan- 
nel and the ionization yield peaks at the zeros of the IR 
field. In Floquet picture, the 13th and 15th HH transition 
matrix elements to the components of 4p Floquet state 
have opposite signs. In general, as the transition matrix 
elements are real, ideally there are only two possibilities; 
either Mq and M2 can have the same sign {(j) = 0, ±27r, ..) 
or opposite signs {(j) = ±7r, ±37r, ..) and we see both vari- 
eties in Fig. ^c). At very low intensity, the 5p ionization 
starts dominating and oscillates with same phase as 2p 
(i.e. 0,±27r). Thus, in our experiment, as the intensity 
decreases or time-delay get longer, we observe a change in 
phase of ion-yield oscillation, going from zero when 2p is 
dominant, through tt where 4p contributes substantially, 
and eventually towards 27r where 5p dominates. 

Our work represents a direct measurement of quan- 
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